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Abstract We have investigated both theoreticahy and 
experimentahy dipolar relaxation in a gas of magnet- 
ically trapped chromium atoms. We have found that 
the large magnetic moment of 6 jub results in an event 
rate coefficient for dipolar relaxation processes of up 
to 3.2 • 10~^^cm^s~^ at a magnetic field of 44 G. We 
present a theoretical model based on pure dipolar cou- 
pling, which predicts dipolar relaxation rates in agree- 
ment with our experimental observations. This very gen- 
eral approach can be applied to a large variety of dipolar 
gases. 



moment of 6/iB heating and atom loss due to dipolar 
relaxation is so significant that it is impossible to reach 
Bose-Einstein condensation by the standard procedure 
of rf-e vapor at ion in a magnetic trap. 

Using a simple model which only takes into account 
the dipolar interaction and neglects all other aspects of 
the molecular interaction on the scattering process we 
are able to explain the observed dipolar relaxation rate in 
ultra-cold ^^Cr. Qualitatively, we find the same behavior 
for ^^Cr. This gives strong evidence that, in the absence 
of resonant behavior, inelastic collisions due to dipolar 
relaxation can be treated in a very general way. 



1 Introduction 

Recently, the dipole-dipole interaction in Bose-Einstein 
condensates has generated significant theoretical inter- 
est It has been shown that the long-range character 
and the anisotropic nature of the interaction between 
magnetic or electric dipoles will lead to many interesting 
new phenomena like novel quantum phase transitions ^ , 
dipolar order pj and spin tunnelling ^ in the conden- 
sate. New questions concerning stability and tun- 
ability of the interaction 7 of such a condensate arise. In 
most of the recent work, the interaction between dipoles 
is treated as completely elastic. However, in a typical ex- 
periment aimed at the creation of a Bose-Einstein con- 
densate, the atoms are held in a magnetic trap in the 
energetically highest magnetic sub-level. Therefore, in- 
elastic dipolar relaxation processes in which the atomic 
spin flips and the Zeeman energy is released cannot a 
priori be neglected. 

Most of the above proposals require a large dipole 
moment like the one of a polar molecule in an electric 
field. In this paper, we show theoretically and exper- 
imentally that already for chromium with a magnetic 



2 Dipolar Relaxation: Theory 

As a consequence of Maxwell's equations it is impossi- 
ble to create a maximum of a static magnetic field in a 
region where external currents are absent. For this rea- 
son atoms are usually magnetically trapped in the en- 
ergetically highest Zeeman sub-level (weak-field-seeking 
state). The trap lifetime of these atoms is mainly lim- 
ited by collisions. Aside from background gas collisions 
and three-body recombination, the atoms undergo spin- 
exchange and dipolar relaxation collisions. In contrast 
to spin-exchange collisions, where the total spin is con- 
served, the spin is flipped in a dipolar relaxation process 
since dipole-dipole interaction does not preserve the to- 
tal spin, but rather the total angular momenta. There- 
fore, if one considers a cloud of fully polarized atoms 
in the energetically highest state, only spin relaxation 
which does not conserve the total spin leads to atom 
loss from the magnetic trap. 

Here, we estimate the inelastic cross sections for spin- 
flip transitions due to the magnetic spin-dipole spin- 
dipole interaction using Fermi's Golden rule. A similar 
form of this approximation, equivalent to the Born ap- 
proximation, has already been considered in ,8^. for 
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calculating the dipolar relaxation rates of meta-stable 
triplet helium. In particular, in £| it is shown that the 
bare Born approximation is reasonably consistent with 
the distorted- wave approximation which is a true first or- 
der calculation. Because of the small dipole-dipole cou- 
pling in polarized triplet helium {S = 1), the latter is 
considered to be an almost exact approach. 

For chromium, the situation is slightly different be- 
cause the magnetic coupling is larger. A first order calcu- 
lation, in the form of the Born approximation, provides 
us with an approximate analytical tool for the interpre- 
tation of the experimental results, but may not be ex- 
pected to be as accurate as in the case of metastable 
helium. 

We consider for simplicity the case of two identical 
atoms trapped in the energetically highest sub-level of 
the Zeeman manifold. Moreover, we consider the case 
of zero electronic angular momenta L = such that 
the magnetic interaction solely arises from the electronic 
spin. The long range part of the interaction between 
two atomic magnetic dipole moments /l^i = gsfJ^B^i and 

= 9sfJ^B^2 located at ri and r2, where Si and S2 are 
the two spin matrices, reads 

, , , ,2 (Si-S2)-3(Si-r) (S2-r) 

^dd(r) = /j.o{gs/j^B) • (1) 

Here we have introduced the interatomic separation r = 
r2— ri with r = r/r and the magnetic permeability of the 
vacuum /io- The tensorial part of the dipolar interaction 
(^, namely (Si • S2) — 3 (Si • r) (S2 • r) can be rewritten 
in terms of spin-flip operators: 



-^(2 i^i, + f_^i+ + f+5i_) 

X (2 zS2z + r_6'2+ + r+6'2-) 



where 



r 



(2) 



(3) 



and where 5'+ = {S^ + iSy) / and S- = {Sx -iSy)/V2 
are ladder operators of the spin operators. From the 
above expression Q it can be seen that in first order 
Born approximation each single atom can flip at most 
one spin {Ams = 1) and in total for the two- atoms sys- 
tem Ams = 2. 

In the Born approximation the total cross section for 
two identical colliding atoms in the same internal state 
can be expressed in the following form 



m 
4^ 



1 

kkf 



|/7dd(k-kO|'^(|k'|-^/)^k^ 



± / /7dd(k-kO/7d*d(-k-kO(^(|k'|-^/)^k' 



(4) 



where k and k' are the initial and final wave vectors, 
respectively, /c/ is the modulus of the final state wave 



vector. t/dd(q) is the Fourier transform of the dipole- 
dipole interaction (^) already contracted between the 
(symmetric) spin (or internal state) functions of initial 
and final states. The sign ± of the second term in the 
square bracket, the "exchange" term, depends on parti- 
cle's symmetry: + holds for bosons and — for fermions. 
When the two colliding atoms are not in the same inter- 
nal state a generalization of with both signs should 
be used. 



2. 1 Dipolar transitions in a dipolar Bose gas 

We have in mind the specific case of ^^Cr atoms in 
the electronic ground state, magnetically trapped in the 
highest Zeeman sub- level. For ^^Cr, both the nuclear 
spin and the electron orbital momentum vanish (1=0, 
L=0), and only the electron spin is finite {S = 3). The 
degeneracy of the 2S+1 atomic ground state sub- levels 
is lifted by a magnetic field B. The Zeeman splitting is 
given by 



= gsfisB 



(5) 



where /i^ is the Bohr magneton and gs is the Lande 
factor {gs = 2). The magnetic field is chosen along the 
z direction and determines the quantization axis. We 
label the internal states via the z-component of the total 
internal angular momentum (in this specific case only 
the electronic spin), like \ms). 

We consider specifically the case in which both atoms 
are prepared in the stretched state, i.e \ms = s). The to- 
tal spin function of the two colliding particles is already 
symmetric and will be denoted as \s]s). The possible 
final internal states of the two colliding atoms are 



1/0) = k; 5) 

\s;s-l) + \s-l;s)) 



1/1) = 
1/2) = 



V2 

s-l;s- 1) 



(6) 
(7) 
(8) 



where the subscripts 0, 1, 2 refer to the number of spins 
that are flipped. For dipolar scattering, the total cross 
section will depend on the relative orientation of the ini- 
tial relative momentum k and the polarization axes. The 
averaged cross sections over all possible orientations for 
the 0, 1 and 2 spin- flip processes reads 



IGtt 4 / /J^o{gsl^B)'^rn^ ^ 



45 



[! + /.(!)] 



(9) 



CTl 



f M9sl^B)^m 



15 



[l^h{kf/h)]f {10) 



While (Ji and (72 are the cross sections for dipolar relax- 
ation, (Jo is the contribution to the elastic cross section 



Dipolar Relaxation in an ultra-cold Gas of magnetically trapped chromium atoms 



3 




Fig. 1 Function li(x); li(x) is a monotonically increasing 
function which asymptotically approaches 1. 



due to the dipole-dipole interaction. Here 



h{x) 



3 jl-x^)^ 
8x(l +^2) 



log 



(12) 



represents the ratio of the exchange contribution to the 
direct one of Eq. Q). The function h (see FigGJ is defined 
in the interval (l,oo) and exhibits a monotonic increase 
from h(l)=-l/2 to h{oo) = 1. (For large values of the 
argument h{x) ^ 1 — A/x'^.) 

Note, that the term in brackets in Eq. (|TH) can be 
expressed in terms of the classical electron radius ro, 
atomic mass m and electronic mass m^: 



l^o{gsl^B)'^rn 



m 



(13) 



One can see that the explicit dependence on ?i disap- 
pears. 

The factor ^ accounts for the different density of 
final states in the inelastic process {AE ^ 0). The ra- 
tio between final and initial wave- vectors is kf/ki = 

^2^2^ for a one spin- flip transition and kf/ki = 



1 



^ for a two spin-flip transition. If the Zeeman 

energy is larger than the initial kinetic energy it exhibits 



a characteristic v B dependence. The expressions in Eq. 
(fTT|) for one and two spin flip processes can also be ap- 
plied to meta-stable helium {S = 1). For B = we 
recover the result of ^^9^ for one and two spin-flip pro- 
cess (using {vrei) = ^l^ksT /irm ). Our method yields 
an analytical solution for all values of the magnetic field 
which reproduces the magnetic-field dependence found 
in IHUni at very small and large magnetic field. 



3 Experiments 

3.1 Dipolar relaxation in chromium 

The magnetically trappable substates in the ^S's ground 
state of chromium are the low- field seeking states ms = 



1,2,3. In a polarized sample of atoms in the energetically 
highest state {ms = 3), there are two possible decay 
channels due to dipolar relaxation: 



|3;3),|3;3) 
|3;3),|3;3) 



|3;3),|3;2) 
|3;2),|3;2) 



AE 
2AE. 



(14) 



In these processes either one or both colliding atoms 
change their state to the ms = 2 state and release a 
Zeeman-energy of AE or 2AE^ respectively. The en- 
ergy is equally distributed among the colliding atoms. 
In a single spin flip transition the increase in tempera- 
ture which corresponds to the released Zeeman energy 
{SksAT = AE/2) is 250 /iK for an offset field of 
5o=10G. This is generally not high enough to eject the 
atoms from the trap. As a consequence thermalization 
leads to heating of the trapped cloud (see Fig.I2^a)). The 
increasing population of the ms = 2 state opens addi- 
tional collision channels due to spin exchange collisions 
and dipolar relaxation between identical and different 
states. The sample becomes more and more depolarized 
and only collisions with a final state of ms < 1 lead to 
atom loss. In this paper, we concentrate on dipolar relax- 
ation of ms = 3 state atoms in an almost spin-polarized 
sample. 



3.2 Methods to determine the dipolar relaxation rate 

In the following, we present three experimental methods 
from which we obtain the dipolar relaxation coefficient 
Pdr for the ms = 3 state. It can be determined by mea- 
suring either the decreasing population (method (i), (ii)) 
of this state or the heating of the atom cloud in the trap 
(method (iii)) over time. Methods (i) and (iii) have also 
been used in JUj to determine the dipolar relaxation rate 
of Cs in the F=3 ground state. 

The atom loss from the ms = 3 state of an initially 
polarized sample of Ns atoms is described by 



dt 



= - J r{r,t)ns{r,t) d^r, 



(15) 



where ns{T^t) and /^(r, t) are the position and time- 
dependent atom density and loss rate, respectively, 
/^(r, t) can be written as 



r(r, t) = jbg + f^drTisir, t) + r(r, t). 



(16) 



Here background gas collisions with a rate coefficient 
75^ are represented by the first term and the density- 
dependent dipolar relaxation with a rate constant Pdr by 
the second term. Other collisional losses e.g. three-body 
recombination, relaxation and spin-exchange processes 
which involve states other than the = 3 state are in- 
cluded in /^(r, t). For low densities and small populations 
in these other states, this term can be neglected. 

The first method (i) takes advantage of the Zeeman- 
energy gain during a relaxation process. The energetic 
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Fig. 2 a) Dipolar relaxation in a magnetic trap. Atoms ini- 
tially prepared in the ms — 3 state undergo dipolar relax- 
ation and fill the ms = 2 state. Indicated is the energy gain of 
a colliding atom of Z\E/2 and Z\E of single and double spin 
flip transition, respectively. Thermalization of the released 
Zeeman energy leads to heating in the trap, b) If a rf-shield 
is applied, the collisional products leave the trap. Thus, the 
collision process is turned into a real loss process. 



collisional products are selectively removed by a radio- 
frequency (rf) shield with frequency Vrf as depicted in 
Fig. 121 Thus, dipolar relaxation is turned into a real loss 
process and can directly be measured by monitoring the 
remaining atoms. The decay is fitted to the solution of 
the differential equation: 



dN^(t) 
dt 



dr 



v{t) ■ 



(17) 



Where V = V8{2Tr)^/'^a^aya^ is the mean volume of an 
atomic cloud with a Gaussian density distribution with 
widths cTx^y^z' To remove the collisional products 
the rf-frequency of the shield has to be chosen such that 
mj{hiyrf —AE) < AE/2. In this case /3'^^ is an event rate 
which contains single and double spin flip transitions. 
P'^^ can be calculated from the thermal average over both 
inelastic collision channels: 



(18) 



The factor of 2 accounts for two removed atoms per 
collision event. At the same time, the rf-frequency has 
to be adjusted to a sufficiently high value, such that 
the initial ms = 3 cloud is not significantly evaporated 
{hurf ^ AE + ksT^ where T is the initial temperature 
of the cloud). From the above considerations it becomes 
clear that for a given temperature this method is only 
applicable for offset fields Bq SksT/fiB' 

For smaller magnetic offset fields, we use a second 
method (ii). The initially polarized cloud evolves in the 
magnetic trap for a certain time t. By applying a strong 
magnetic field gradient after switching off the trap, the 
magnetically trapped atoms can be spatially separated 
with respect to their substates before imaging. This is 
equivalent to a Stern- Gerlach experiment. A series of 
these images yields the time evolution of the population 



of the trappable sub- levels. From the initial decay of the 
ms = 3 we obtain the loss rate Pdr- 



dNsjt) 
dt 



V{t) ' 



(19) 



It is not possible to distinguish between single (one atom 
lost from ms = 3) and double (two atoms lost from 
ms — 3) spin flip transitions with this method. Hence 
the relation between and the cross sections changes 
to 

/3,, = ((ai+2-a2)^.e/).._. (20) 

The practicability of this method is limited by the mag- 
nitude of the magnetic field gradient which can be ap- 
plied to separate the expanding thermal clouds before 
imaging. 

The third method (iii) measures the heating of the 
atomic cloud due to the energy released by dipolar re- 
laxation processes. The evolution of the temperature is 
given by 



~ Pdr TT/ • 



dt 



V{t) 



(21) 



Since a two spin flip transition yields twice the energy 
of a single spin flip transition the obtained value is Pdr 
as in method (ii) fEa l20j) . This method requires that the 
cloud is in thermal equilibrium, i.e. that rethermalization 
by elastic collisions is faster than the dipolar relaxation 
process itself. As the velocity of atoms which have un- 
dergone a dipolar relaxation event and the volume which 
is occupied by the collisional products increases with in- 
creasing offset field, the elastic collision rate drops and 
impairs the measured result. 



3. 3 Preparation of the sample and experimental 
methods 

The basic setup used in our experiments has been de- 
scribed previously PTUT^ and is only briefly summa- 
rized in the following. Zeeman-slowed ^^Cr atoms are 
magneto-optically trapped in the magnetic field configu- 
ration of a loffe-Pritchard trap in the cloverleaf configu- 
ration using the ^S3-^P4 transition. Since this transition 
is not closed, we populate the metastable state via 
the weak ^P4-^D4 intercombination line. In this state 
the atoms are decoupled from the trapping light and 
can be magnetically trapped because of the high mag- 
netic moment (Q/j^b)- We load about 2-10^ atoms con- 
tinuously in the loffe-Pritchard type trap ("CLIP" -trap 
jl3j). After switching off the trapping light, we opti- 
cally transfer (using light at 658 nm) the atoms back 
to the ^Ss ground state, which possesses the same mag- 
netic moment (/i = Q/j^b)- Subsequently, we compress the 
sample and perform Doppler cooling in the magnetic 
trap. Using this method, we obtain a nearly polarized 
{ms = 3) atomic ensemble with an initial density of 
about 10^^ atoms/cm^ and a temperature of 275 ±25 jaK. 
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This is the starting point for most of our measure- 
ments of dipolar relaxation at high magnetic offset fields 
(17 — 44 G) using method (i). However, for smaller off- 
set fields (<25G), it was necessary to reduce the tem- 
perature further by rf-evaporation to realize a cut-off 
parameter for the rf-shield of 77 > 5 14 . After prepara- 
tion, the rf-shield is applied and the atomic ensemble is 
allowed to evolve in the trap for a variable time t. The rf- 
frequency is set to 1 MHz below 7AE/6h. Subsequently, 
the shadow of the atoms is imaged onto a CCD-camera 
by illuminating the atoms for 100 jas with a circularly 
polarized beam after 5 ms of free ballistic expansion in a 
homogenous magnetic field. 

The number of remaining atoms Ns{t) and the mean 
volume V{t) are extracted from a 2D-Gauss fit to ab- 
sorption images using the separately measured trap fre- 
quencies (at 27 G offset field: Vradiai = 120 Hz, Uaxiai = 
73 Hz). A series of these measurements with different 
evolution times t yields Ns{t). The temperature T{t) is 
obtained by recording the ballistic expansion of the atom 
cloud at different evolution times. 

For measurements at low offset fields (<2 G), we use 
the second approach (ii). Since the temperature obtained 
after Doppler cooling is too high to separate the atoms 
with the available field gradient, we reduce the tempera- 
ture of the sample to 50 /iK by rf-evaporative cooling at 
an offset field of 0.7 G (trap frequencies for this config- 
uration are typically ^radial = 806 Hz, Uaxiai = 42 Hz). 
Subsequently, we turn off the rf-frequency and change 
the offset field for a free evolution time t to the de- 
sired value. Shortly (10 ms) before the Stern-Gerlach se- 
quence, the field is changed back to 0.7 G. 

The separation of the substates is accomplished by 
switching off the axially confining pinch and bias fields 
of our clover-leaf trap and simultaneously applying a ho- 
mogenous field in the imaging direction which is perpen- 
dicular to the offset field direction and diagonal to the 
leaf position of our trap. By this means the center of 
the trap is displaced radially and the atoms get acceler- 
ated in the gradient field of the leaves (170G/cm) along 
an axis perpendicular to the imaging direction. After a 
1.7 ms acceleration phase, all magnetic fields are turned 
off and the atoms expand freely for another 5 ms before 
the spatially separated clouds are imaged on the CCD 
camera by illuminating them for 1 ms with a calibrated 
light field. The data analysis is perform by averaging 
three pictures taken under identical experimental condi- 
tions for each evolution time. An example for such an 
image taken at 0.3 G is depicted in Fig. 01 The clouds 
have been accelerated downwards. The cloud containing 
atoms in the energetically highest state {ms = 3) there- 
fore appears at the lowest position. 

The number of atoms in each substate are obtained 
by 2D-Gauss fits within an equally sized rectangular area 
around each cloud. Because of the long exposure time 
and the trap geometry being changed while switching 
off the trap, temperature and volume were obtained in 



m3=1 
ms=2 

E 

:: 

Fig. 3 Typical fluorescence picture taken after separation of 
the sub- levels in a Stern-Gerlach type experiment and 5 ms 
time of flight. The sub- levels have been accelerated down- 
wards. 

separate measurements where we record the free expan- 
sion of the whole cloud at different evolution times with 
a exposure time of 200 /is. 

We determine the offset field by probing the min- 
imum of the trapping field via rf removal of trapped 
atoms. The achieved accuracy of the values is about 
300 mC for measurements at high offset values and 
100 mC for the measurements at low offset fields. The 
background gas collision rate 75^ = 1/200 s is obtained 
from a separate measurement with a low density atomic 
sample. This reduces the fitting parameters to the initial 
number of atoms and /3dr- The increasing volume due to 
heating during the measurement was included in the fit. 
The errors in the determination of the number of atoms 
is about 30%. The error bars shown in the figures are 
the square root of the diagonal elements of the covari- 
ance matrix for the fitting parameters obtained from a 
least-squares fit. 

3.4 Results 

We have performed measurements using the first method 
(i) for different magnetic offset fields ranging from 17 to 
52 G. An example of a typical decay curve of the popu- 
lation Ns{t) of the ms = 3 level, taken at an offset field 
of 27 G, is plotted in Fig. El The fit to a two body decay 
using the solution of Eq. (|T7)) yields Pdr = (2.5 ± 0.5) • 
10^-^ cm^/s and is indicated by the black curve in the fig- 
ure. The inset to this figure demonstrates the evolution 
of the temperature with (bullets) and without (squares) 
rf-shield. Hereby, the temperature is derived from the 
width of the cloud. Without rf-shield, strong heating 
with an initial rate of more than 600/iK/s is observed. 
With rf-shield, the temperature remains constant. 

The results of the decay measurements are summa- 
rized in Fig. [51 The scatter of the data for a given offset 
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Fig. 4 Decay measurement using method (i). The number 
of atoms A^3(t) remaining in the ms = 3 state after a vari- 
able free evolution time in the trap at an offset field of 27 G 
is depicted. The curve represents a fit using the solution of 
Eq. El The inset provides the temperature evolution with 
(bullets) and without (squares) rf-shield. 
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Fig. 5 Dependence of the dipolar rate constant p^r on the 
offset field Bq. Depicted are p^r (bullets), the calculated val- 
ues (diamonds, connected by straight lines to guide the eye) 
and values obtain from a different setup (open squares) . The 
calculated values are obtained by using Eq. (UHl and the ex- 
perimental parameters T and Bq 



field indicates the uncertainty in our measurement. As 
expected, the loss coefficient increases with increasing 
field. Each bullet represents a separate decay measure- 
ment. Also included are two values (depicted by open 
squares), which were measured in a different experimen- 
tal setup using the same method. The diamonds which 
are linked by a straight lines to guide the eye, represent 
the results obtained from the theoretical model using 
Eq. (|18|) and the experimental parameters T and ^o- 
Though there is a slight discrepancy between the data 
and the model calculation, the experimental values are 
reproduced within a factor of two. The deviation be- 
tween the two experimental data sets is attributed to 
systematic errors in the determination of the density in 
these experiments. 

Measurements of the dipolar relaxation rate at low 
offset fields have been performed using method (ii). Fig. 
ini shows a typical result of a Stern- Gerlach Experiment 
at an offset field of 0.7 G. The figure contains the evo- 
lution of the number of atoms in each sub-level and the 
total number of atoms {Ntot{t) = Ni{t) ^ N2{t) ^ Ns{t)) 
during the first 15 s. While the total number of trapped 
atoms (circles) stays approximately constant, about 35% 
of the atoms initially in the ms = 3 state (squares) get 
redistributed among the other trapped states. The ini- 
tially prepared atomic cloud is not fully polarized even at 
t = s. We observe a small fraction of the cloud (about 
20%) in the ms = 1,2 states (triangles, diamonds). The 
origin of this population remains unclear and might ei- 
ther be caused by incomplete rf evaporation or by Ma- 
jor ana spin flips during switching of the magnetic 
fields. An indication that evaporative cooling populates 
these states is the slightly decreasing number of atoms 



x10' 




5 10 

time [s] 



15 



Fig. 6 Redistribution among the sub-levels. While the total 
number of atoms (A^tot) (circles) stays constant, the popula- 
tion of the initially almost polarized sample gets redistributed 
by dipolar relaxation and spin exchange collisions. The pop- 
ulation NmS of each sub- levels ms = 1,2 and 3 is indicated 
by triangles, diamonds and squares, respectively. Also shown 
are the fitting curves to the decay of the ms = 3 (dotted 
line) and loading of the ms = 2 state (dashed line). 



which can be observed in the first few hundred ms in the 
ms = 1,2 states in samples with lower densities. 

The decay of the ms = 3 population is non-exponen- 
tial. The dynamics is described by two coupled, non- 
linear differential Eqs. (ITUj) and (PT|) (where V{t) ex 
T^/^(t)). By assuming a linear increase in the volume 
which is describing by our data well, we obtain a solution 
for Eq. (|T^ . A fit to data for which the population of the 
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offset field [G] 

Fig. 7 Dependence of the dipolar relaxation rate [3dr on the 
offset field Bq at low magnetic fields. The experimental val- 
ues obtained by method (ii) are indicated by bullets. The 
diamonds connected by straight lines to guide the eye are 
calculated using Eq. and experimental parameters Bq 
and T. Also included are the values (triangles) obtained from 
the temperature increase of the sample. 



ms = 2 state is below 25% yields a dipolar relaxation 
rate constant oi pdr = (3.1 ± 0.4) • IQ-^^cm^s-^ The 
fit to N3(t) is illustrated by a dotted curve. We repeated 
this measurement for 5 different offset fields ranging from 
0.15 to 1.7 G. The dependence of the measured values on 
the offset field is shown in Fig. [7| Again, the diamonds 
which are connected by straight lines to guide the eye 
are obtained from Eq. using the measured experi- 
mental parameters Bq and T with no other adjustable 
parameter. The agreement of the calculated values with 
the experimental data is reasonable, although the square 
root dependence of the decay rate on the magnetic field 
is not exactly reproduced. 

We have also investigated the loading process of the 
ms = 2 state in order to deduce the two body-loss coeffi- 
cient (32 of this state. For short lifetimes, the loading rate 
into the ms = 2 state can be exclusively attributed to 
the loss out of the ms = 3 state. Neglecting background 
gas collisions, we obtain the following rate equation: 



dN2{t) dNs{t) 



dt 



dt 



•A 



V2{t) 



(22) 



where V2 = (3/2) ^/^Vs is the mean volume of the ther- 
malized ms = 2 cloud. As spin exchange collisions and 
dipolar relaxation are possible in this state, P2 contains 
both processes. Assuming a constant volume, the equa- 
tion can be solved. A fit of the first order power series 
expansion of the solution to our data yields a value of 
(32 = (1.1 ± 0.2) • 10~^^ cm^5~^, where we have taken the 
average of the measurements between 0.1 and 1.7 G. The 
fit is indicated in Fig.EJby a dashed curve. The obtained 
values are shown in the inset of Fig.[3 The values of the 




6 8 10 
time [s] 

Fig. 8 Heating due to dipolar relaxation at a low offset field 
(0.7 G). The solid line is a linear fit to the data. 



coefficients are about a factor 50 higher than the dipo- 
lar relaxation coefficient of the ms = 3 state. This can 
probably be attributed to the additional decay channels 
mentioned in SecEl 

We observe that the population in the ms = 1 state 
increases initially much slower than the population in the 
ms = 2 state. This indicates that ms = 3 state does not 
directly decay to the ms = I state via dipoar relaxation. 
The almost constant total number of atoms gives further 
evidence that the observed loss from the ms = 3 state 
involves mainly spin- flips with Ams = 1. Thus, the first 
order Born approximation is a valid approach to describe 
dipolar relaxation in our system. 

An example for the time evolution of the mean tem- 
perature is depicted in Fig. [SI Even at a rather low off- 
set field of 0.7 G the temperature of initially at 50 /iK 
increases rapidly with an approximate rate of 3/iKs~^. 
We estimate the dipolar relaxation rate constant using 
method (iii) by solving Eq. (|2U using the solution for the 
number of atoms in the ms = 3 state obtained above. 
The results of a fit to the linearized solution are included 
in Fig. 13 (triangles) and are in agreement with the results 
using method (ii). 

The dipolar relaxation rate is calculated supposing 
that the inelastic process is dominated by the dipole- 
dipole interaction. Details of the molecular interaction 
are not taken into account. We have tested this assump- 
tion by comparing the heating of the atomic cloud due 
to dipolar relaxation for the two isotopes ^^Cr and ^^Cr. 
The temperature evolution of the atomic clouds at an 
offset field of Bq = 20 G is shown in Fig. El Here, the 
temperature was derived from the 1 / width ay of the 
cloud. Equally high heating rates (13/iKs~^) yielding 
the same loss coefficient can be observed for both clouds 
prepared with the same initial atom number density of 
10^^ atoms/cm^. The data shown in Fig. El indicate that 
the loss due to dipolar relaxation is independent of the 
details of the interatomic potential. Since the scattering 
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time [s] 

Fig. 9 Heating rate due to dipolar relaxation at 20 G off- 
set field for the two isotopes ^^Cr (open circles) and ^°Cr 
(squares). The initial density for both isotopes is roughly 
no = 10^° cm-^ 

length for a{^^Cr) = 50 ao and a{^'^Cr) = 170 ao dif- 
fers significantly 12 one would expect that any effect 
due to the 'molecular' interaction also differs for the two 
isotopes. 



4 Conclusion 

We have presented three different methods to obtain the 
dipolar relaxation rate constant in magnetically trapped 
clouds of chromium atoms. Using these methods we have 
measured the dependence of the rate constant on the 
magnetic offset field of the trapping potential. We find a 
typical dipolar relaxation rate constants of /3dr = 4 • 
10~^^cm^s~^ at a magnetic offset field of Bq = 1 G. 
Thus, even at low magnetic fields, the loss coefficient 
for chromium is one order of magnitude larger than the 
one reported for Gs in lowest hyperfine level (F=3) jlOj 
and another 3 orders of magnitude larger than in Na 

A theoretical estimation of the inelastic cross sections 
for spin-flip transitions due to the magnetic dipole-dipole 
interaction, predicts cross sections for single and double 
spin flip transition which scale with the total angular mo- 
mentum to the third and second power, respectively. The 
values predicted by the calculation are in good agree- 
ment with the experimentally observed values. The influ- 
ence of the molecular potential on the dipolar relaxation 
cross section for chromium was found to be insignificant. 

Our findings could be especially important for ex- 
periments which aim at producing high density samples 
of atoms and molecules with high magnetic or electric 
dipole moments in magnetostatic or electrostatic traps 
pTHTSlT lQ . Large inelastic two-body collision rates are to 
be expected with increasing dipole-dipole interaction, if 
atoms are not trapped in the lowest energy state of the 
dipole. 



As has already been pointed out in ^Hl dipolar relax- 
ation can be almost totally suppressed by polarizing the 
sample in the energetic lowest spin state, which requires 
different trapping schemes. In this state neither dipo- 
lar relaxation nor spin exchange collisions are possible 
at reasonable offset field. The only expected loss mecha- 
nism in this case is three-body recombination. However, 
because of the long-range nature of the dipole potential 
the latter could also be increased compared to atoms 
with a small dipole moment. As a quite general conclu- 
sion, strong dipolar gases can only be cooled to degen- 
eracy in traps that keep the dipole moment in its lowest 
energy state. 

Our results, though not promising at first glance, will 
help us to develop a strategy for reaching Bose-Einstein 
condensation even under the special conditions set by the 
large magnetic moment of chromium. In cesium, losses 
due to inelastic collisions in a magnetic trap have pre- 
vented Bose-Einstein condensation for a long time in a 
similar way as in the chromium case. However, a care- 
ful study of the nature of the inelastic losses has finally 
allowed to devise a successful strategy for condensation 

m- 

This work was supported by the Deutsche Forschungs- 
gemeinschaft and the RTN network 'Gold Quantum 
Gases'. 



References 

1. M. Baranov, L. Dobrek, K. Goral, L. Santos, M. Lewen- 
stein: Phys. Scripta T102, 74 (2002)and reference therein 

2. K. Goral, L. Santos, M. Lewenstein: Phys. Rev. Lett. 88, 
170406 (2002) 

3. H. Pu, W. Zhang, P. Meystre: Phys. Rev. Lett. 87, 140405 
(2001) 

4. H. Pu, W. Zhang, P. Meystre: Phys. Rev. Lett. 87, 
0904001(2002) 

5. K. Goral, K. Rzazewski, T. Pfau: Phys. Rev. A 61, 051601 
(2000) 

6. S. Yi, L. You: Phys. Rev. A 63, 053607 (2001) 

7. S. Giovanazzi, A. Gorlitz, T. Pfau: Phys. Rev Lett. 89, 
130401 (2002) 

8. G.V. Shlyapnikov, J.T.M. Walraven, U.M.Rahmanov, 
M.W. Reynolds: Phys. Rev. Lett. 73, 3247 (1994). 

9. P.O.Fedichev, M.W. Reynolds, U.M. Rahmanov, G.V. 
Shlyapnikov: Phys. Rev. A 53, 1447 (1996). 

10. D. Guery-Odelin, J. Soding, P. Desbiolles, J. Dalibard: 
Europhys. Lett. 44(1), 25-30 (1998) 

11. P.O. Schmidt, S. Hensler, J.Werner, T. Binhammer, A. 
Gorlitz, T. Pfau: J. Opt. Soc. Am. B 20(5), 960 (2003) 

12. P.O. Schmidt, S. Hensler, J.Werner, A. Griesmaier, A. 
Gorlitz, T. Pfau, A. Simoni: |ar Xiv:qua nt-ph/0303069 v2 

13. P.O. Schmidt, S. Hensler, J.Werner, T. Binhammer, A. 
Gorlitz, T. Pfau: J. Opt. B 5, S170 (2003) 

14. N. Masuhara, J. Dolye, J.G. Sandberg, D. Kleppner, T.J. 
Greytak, H.F. Hess, G.P. Kochanski: Phys. Rev. Lett. 61, 
935 (1988) 

15. W. Petrich, M.H. Anderson, J.R. Ensher, E.A. Gornell: 
Phys. Rev. Lett. 74, 3352 (1995) 



Dipolar Relaxation in an ultra-cold Gas of magnetically trapped chromium atoms 



9 



16. A. Gorlitz, T.L. Gistavson, A.E. Leanhardt, R. Low, A. P. 
Chikkatur, S. Gupta, S. Inouye, D.E Pritchard, W. Ket- 
terle: Phys. Rev. Lett. 90, 090401 (2003) 

17. J.D. Weinstein, R. deCarvalho, J. Kim, D. Patterson, B. 
Friedrich, J.M. Doyle: Phys. Rev. A 57, R3173 (1998) 

18. J.D. Weinstein, R. deCarvalho, T. Guillet, B. Friedrich, 
J.M. Doyle: Nature 395, 148 (1998) 

19. H.L. Bethlem, G. Berden, F.M.H Crompvoets, R.T. 
Jongma, A. J. A. van Rolj, G. Meijer: Nature 406, 491 
(2000) 

20. T. Weber, J. Herbig, M. Mark, H.-C. Nagerl, R. Grimm: 
Science 299, 232 (2003) 



